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Abstract 



We study a dilute Bose gas of atoms whose scattering length o is large 
compared to the range of their interaction. We calculate the energy density £ 
of a homogeneous Bose-Einstein condensate (BEC) to second order in the low- 
density expansion, expressing it in terms of a and a second parameter A* that 
determines the low-energy observables in the 3-body sector. The second-order 
correction to £ has a small imaginary part that reflects the instability due to 
3-body recombination. In the case of a trapped BEC with large negative a, 
we calculate the coefficient of the 3-body mean-field term in £ in terms of a 
and A=K. It can be very large if there is an Efimov state near threshold. 



Bose-Einstein condensates of atoms have been extensively studied both experimentally 
and theoretically for 2 cases: liquid helium |l|] and dilute vapors of alkali atoms The 
low-energy length scale i = {mCe/h'^Y^'^ set by the van der Waals interaction — Cg/r® is the 
natural length scale for the S-wave scattering length a and other parameters in the low-energy 
expansions of scattering amplitudes. In the case of liquid helium, the interparticle spacing 
n~^^^ = 3.6 A is much smaller than £ = 37 A and a = 104 A 0. Because the diluteness 
variable y/nP ~ 33 is much greater than 1, the properties of the Bose-Einstein condensate 
(BEC) in liquid helium depend on the detailed behavior of the interatomic potential. In 
contrast, BEC's consisting of dilute gases of atoms with ni^ <C 1 have universal properties 
that depend on the interatomic potential only through the single low-energy parameter a. 
For example, the first few terms in the low-density expansions for the energy density and 
the condensate fraction of a homogeneous condensate can be calculated as expansions in 
powers of \fna?. 

A fundamental open problem is the behavior of a dilute BEC {n£^ <^ 1) with large 
scattering length (a ^ i) when na^ is comparable to or much greater than 1. The diluteness 
condition ni^ <^ 1 excludes the case of liquid helium. The condition na^ ~ 1 implies 
that the low-density expansion is not applicable, because the expansion parameter Vi67ma^ 
is not small. The most basic question is whether a BEC even exists as a well-defined 
quasistable state in this limit. If so, does it have any universal properties that depend on 
the interatomic potential only through the scattering length a? These questions can be 
investigated experimentally in vapors of alkali atoms by tuning a background magnetic field 
to a Feshbach resonance |Q. They can also be studied using numerical methods 

In this letter, we take a small step towards addressing this problem by studying the 
homogeneous BEC with large scattering length (a ^ i) in the extremely dilute limit na^ -C 
1. We calculate its energy density to second order in Vna^, determining the coefficient of 
na^ in terms of a low-energy parameter A*. The parameter A* can be determined from any 
low-energy 3-body observable, e.g. the binding energy of the shallowest 3-body bound state. 
We also determine the coefficient of the 3-body term in the mean-field contribution to the 
energy density for a < 0. Finally we comment on the implications of our results for the case 
of large na^. 

The dilute Bose gas with large scattering length is distinguished from the generic case 
by having a 2-body scattering amplitude / = —a/{l + iak) that has the scale- invariant form 
f ^ i/k for wavenumbers in the range l/a-cA;<^l/£|^. Ifa>0, itis also characterized 
by the existence of a 2-body bound state with binding energy B2 = h'^/ma'^. As first pointed 
out by Efimov |Q , the 3-body sector also exhibits universal properties in the limit of large 
scattering length, such as the existence of a large number of 3-body bound states. The 
number of these Efimov states is roughly \n{\a\/i)/iT. The Efimov spectrum and other low- 
energy 3-body observables involving wavenumbers k <^ 1/i are universal in the sense that 
they depend on the interatomic potential only through a and a single 3-body parameter. A 
convenient choice for this 3-body parameter is the low-energy parameter A* introduced in 
Ref. i. 

In the dilute limit na^ <^ 1, the properties of the Bose gas can be calculated using 
the low-density expansion. The first few terms in the low-density expansion of the energy 
density of a homogeneous BEC are 
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£ = 2-nh'^an'^/m[l + 128/(15 
+ 8(47r - 3^3) /3 [in(na^) + 4.72 + 2B\ na^ + ...}. (1) 



The yno? correction was first calculated by Lee and Yang in 1957 |10|. The coefficient of 



the logarithm in the no? correction was calculated in 1959 [|rT|. The constant under the 



logarithm was calculated by Braaten and Nieto in 1999 ||T2[. It was expressed in terms 
of an effective 3-body coupling constant gzi^) that can be determined by measuring the 
low-energy behavior of the 3-atom elastic scattering rate. This "running coupling constant" 
depends on an arbitrary wavenumber k: 

gsin) = 3847r(47r - 3\/3)[ln(Ka) + B]h^a^/m, (2) 

where B is the same constant as in (|l|). The na^ correction in (Q) is the sum of a mean- field 
contribution g3{K,)n^ /36 coming from the effective 3-body contact interaction and a 2-loop 
contribution from quantum fluctuations around the mean field, which also depends on k. The 
dependence on k cancels in (|l|), reflecting the arbitrariness in the separation of the energy 
density into mean-field and quantum-fluctuation contributions. The next nonuniversal term 



in the low-density expansion is determined by the effective range for 2-body scattering []T3 
It is suppressed by (na^^^'^i/a, and is therefore completely negligible in the limit of large 
scattering length. 

To determine the constant B in the expression (|ID for the energy density for an extremely 
dilute Bose gas with large scattering length, we need to calculate gsi^n) as a function of a 
and the low-energy 3-body parameter A,, introduced in Ref. This can be accomplished 
by calculating the T- matrix element for 3-atom elastic scattering in the low-energy limit. 
As the total energy E of the 3 atoms goes to 0, the T-matrix element is the sum of divergent 
terms proportional to 1/E, 1/ a/E, and In E |]l4l and a remainder. The dependence of the 



T-matrix element on A* enters only through the remainder, which includes a term —g^^n). 
For the case a < 0, Efimov used simple probability arguments to deduce the dependence of 
the remainder on A* [Q. These arguments imply 

5 = 6i + 62tan(soln(|a|A,)+/5) (a < 0), (3) 

where Sq = 1.0064. Efimov did not determine the functional form of B for a > 0. However 
unitarity requires that it have an imaginary part that is related to the 3-body recombination 
rate into the shallow 2-body bound state [|T^,|r^ : 

Im5 = 0.022cos2(soln(aA,) + 1.76) (a > 0). (4) 

The 3-body recombination rate has zeroes at values of aA^, that differ by multiplicative 
factors of 22.7, and consequently Imi? also vanishes at those points. 

A convenient way to calculate 5'3(k), both for a > and a < 0, is to use the effective field 
theory method developed by Bedaque, Hammer, and van Kolck in Ref. 0. The effective 
theory is defined by the lagrangian density 
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where we have set ^ = 1. The auxihary field d, which annihilates a pair of atoms at a point, 
can be eliminated using the equation of motion d = + mg^/ {727ra)\ip\'^ + ...). However, 
it is more convenient to keep it, because 3-body observables can be conveniently calculated 
in terms of 4-point Green's functions of i/), d, ijj^ and d^. For example, the T-matrix element 
for 3-atom elastic scattering can be expressed as a sum of 9 connected 4-point Green's 
functions. If the incoming atoms have momenta ki, k2, and ks with ki + k2 + ka = 0, the 
external energies and momenta of ip and d are set to (fc^/2m, ki) and ((A;2 + fc|)/2m, k2-|-k3), 
respectively, and similarly for the outgoing atoms. The Green's function is then summed 
over the 3 cyclic permutations of ki, k2, and ks, and over the 3 cyclic permutations of 
the final momenta to get the T-matrix element. The Green's functions can be calculated 
nonperturbatively by solving the integral equations shown in Fig. la. The single lines are 
propagators for ip and the double lines are exact propagators for d. An ultraviolet cutoff 
A must be imposed on the loop momentum in the integral equation. The 3-body coupling 
constant in is tuned as a function of A so that the low-energy 3-body observables 
depend only on a and the low-energy parameter A^. 
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FIG. 1. (a) Integral equation satisfied by the off-shell amplitude T{p). (b) Definition of the 
subtracted amplitude T{p). 

It is not necessary to calculate a physical observable to determine (73 (k). For small 
external energies and momenta {00, k) with m|ci;|a^ -C 1 and |k| -C 1, a Green's function 
associated with the lagrangian @ can also be calculated perturbatively in a and g^ using the 



methods of Ref. |]T2|. By computing the same Green's function using the nonperturbative 
method of Ref. we can determine g^ as a function of a and A^.. We choose the two particle- 
irreducible connected Green's function T{p) with external energies and momenta (0, 0) for 
and d, (p^/2m, p) for tp^ and (— p^/2m, — p) for d"''. As p ^ 0, T{p) has divergent terms 
proportional to 1/p, and ln{p) that come from the diagrams in brackets in Fig. lb. 

The coupling constant (73 (k) in Ref. [0 was defined using dimensional regularization in 
3 — 2e dimensions to regularize the ultraviolet divergences and minimal subtraction with 
renormalization scale k to remove the poles in e. Calculating the diagrams that contribute 
to T{p) in the limit p using perturbation theory to first order g^ and 4th order in a, we 
obtain 

T{p) — > -^3(k)/36 + 16n^aym[l/{p'^a^) - 2%/ (Spa) 

-2(47r - 3v^)/(37r) Mp/k) + Ams] }, (6) 
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where Ams = —0.77268 and the subscript MS stands for "minimal subtraction." The non- 
perturbative calculation of the Green's function T(p) using the method of Ref. ^ requires 
solving the integral equation in Fig. la. It is convenient to define a subtracted Green func- 
tion T{p) that has a finite limit as p — using the subtraction shown in Fig. lb. We solve 
the integral equation for T{p) as a function of p, a, and A,, and then set p = 0. The resulting 
expression for T{p) in the limit p — > is 

T{p) — > I67i^a'^/m{l/{p'^a^) - 2n/{3pa) 

-2(47r - 3V3)/(37r) [\n{pa) + A{aA,)] }, (7) 

where A{aA^) is a function of aA^, that is determined numerically. Matching the expressions 
(^) and (1^), we find that the constant B in (2) is given by -B = —Ams + A{aA^). 

We first consider the case of large positive scattering length. For a > 0, our result for B 
is complex. Its imaginary part is equal to (^), as required by unitarity. For the real part, 
we obtain a good fit with the empirical formula 

Re S = 1.22 + 0.021 cos2(soln(aA,) + 1.0) (a > 0). (8) 

The real part is much larger than the imaginary part given in (p. The oscillatory term in 
the real part has almost the same amplitude as the imaginary part, but a different phase. 

In the coefficient of na^ in the expression (|l|) for the energy density, the real part of 
the factor in brackets can be separated into 2 terms: 21n(Ka) + 2Re5, which comes from 
the mean-field contribution to the energy density, and ln(na/K^) + 4.72, which comes from 
quantum field fiuctuations. Although k is arbitrary, we can by a suitable choice of k arrange 
for most of the na^ correction to come from the mean-field contribution. The smallest value 
of K that is physically reasonable is of order ^/l67^m, which is the inverse of the coherence 
length. Below this scale, the dispersion relation for the atoms is modified by collective effects 
and 3-atom scattering is therefore no longer described accurately by the vacuum T-matrix. 
If we choose k = 1.5\/l67ma, the entire na^ correction in (|l]) is taken into account by the 
mean-field contribution g3{K,)n^/36. 

The na^ correction to the energy density has an imaginary part that comes from the 
imaginary part of B given in (^). This imaginary part refiects the fact that the BEG is 
only a quasistable state. Atoms will be continually lost from the condensate by 3-body 
recombination into the shallow 2-body bound state and a recoiling atom. In the mean field 
approximation, the rate of decrease of the number density is {d/dt)n = —{lmg2/6h)n^. Note 
that Irag^ is independent of k. The energy loss rate is given by dS/dt = {dS /dn)dn/dt. 
The condensate is particularly stable at values of a that correspond to ln(aA*) = 2.93 mod 
71 /sq. At these discrete values of a that differ by multiplicative factors of 22.7, there is 



an interference effect that causes the 3-body recombination rate to vanish |]I5|,|T6|. The 
recombination process that gives the energy density an imaginary part produces atoms and 
bound states with wavenumber k = 2/{^/3a). These high energy atoms and bound states 
can subsequently thermalize and cool by elastic scattering with atoms in the condensate. 
The cross section for elastic scattering of the atoms and bound states can be calculated as a 
function of a and A^,. For example, the scattering length was calculated in Ref. and the 
effective range in Ref. [0. 
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We now turn to the case of large negative scattering length. If a < 0, the expression 
(HI) for the energy density does not apply, because the homogeneous condensate is unstable 
to collapse. However a condensate with a < can be stabilized by a trapping potential as 
long as the number of atoms is below some critical value A^max |18|- An effective 3-body 
contact interaction —g^ gives a mean-field contribution gsn^/ 36 to the energy density. The 
effect of such a term on the stability of a trapped condensate has been studied by adding 
a 3-body term to the Gross-Pitaevski equation [^. Even a small positive value of can 
considerably increase the critical number A^max- 

In the case of large negative scattering length, the 3-body coupling constant (73 (k) is given 
by (|^). The dependence of our result for B on aA* has the functional form predicted 
by Efimov. The values of the constants are bi = 1.23, 62 = —3.16, and (3 = 0.19. The 
appropriate choice for k in this case is the lowest wavenumber above which 3-body scattering 
can be described by the T-matrix for scattering in the vacuum. The vacuum T-matrix 
becomes inaccurate not only because of the nontrivial quasiparticle dispersion relation but 
also because of the inhomogeneities of the trapping potential V{r). The minimum value 
of K set by the coherence length is > Am{fi — V{r)), where n is the chemical potential. 
For a harmonic potential with frequency u, the minimum set by the inhomogeneity of the 
potential is > muj/fi. In previous studies of the mean- field effects of the 3-body term Jiy 



the coefficient g^ was assumed to be a constant. If the large scattering length is obtained 
by tuning the magnetic field to a Feshbach resonance, g^ is given by the expression @, 
which scales roughly like a^. Note that the constant B given in @ diverges at values of 
a that correspond to ln(aA=K) = 1.37 mod tt/sq. At these discrete values of a, there is an 
Efimov state at the 3-atom threshold. Near these values of a, the 3-body term in the Gross- 
Pitaevski equation becomes particularly important. However if the Efimov state is too close 
to threshold, the mean-field approximation breaks down because the energy dependence 
of the 3-body elastic scattering amplitude from terms proportional to 1/{E + B^) is not 
properly taken into account. 

We turn finally to the problem of a homogeneous BEG with a large scattering length 
(a ^ i) that is dilute enough so that ni^ -C 1 but dense enough so that na^ ^ 1. Our 
results for the extremely dilute limit provide a reference point for discussing this problem. 
The fact that our energy density has an imaginary part serves to emphasize that such a 
state will at best be quasistable. Monte Garlo methods that search for the absolute ground 
state will therefore be of limited utility. In the extremely dilute limit, the condensate is 
quasistable because the imaginary part of the energy density is suppressed by na^. Is there 
any mechanism that can provide quasistability at large na^? If a quasistable Bose- Einstein 
condensate does exist, there are some other obvious questions. At what rate does the number 
density in the condensate decay, and what is the fate of the atoms that disappear from the 
condensate? We have answered these questions for the extremely dilute limit. 

Assuming that a quasistable BEG exists in the limit na^ —>■ 00, does it have universal 
properties that depend only on a? If so, dimensional analysis implies that the energy density 
must have the form £ = Cfi^n^^^ /m, where C is a constant. In the extremely dilute limit, 
we found that £ depends on A* at second order in yna?. Does £ depend on A* in the limit 
no? — >• 00? If so, the coefficient C in the energy density is not a constant, but depends on 
n in a peculiar way. It must be a periodic function of ln(n~^/^A^,) that returns to the same 
value when n is increased by a factor of about 11,700. This follows from the discrete scaling 
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symmetry of low-energy 3-body observables discovered by Efimov [|8|,[T7[], which imphes that 
C{n~^/^K^) = C{22.7n~^^^A^). This discrete scahng symmetry is related to the fact that 
there are Efimov states with sizes differing by factors of 22.7 and ranging all the way from 
order i to order a. The question of whether C depends on A^, is not addressed by the 
constrained variational calculations of Ref . . Their variational ansatz excludes states with 

3- body configurations that are sensitive to A^,. For example, in the 3-body sector, it would 
exclude the Efimov states. 

If the coefficient C in the energy density depends on the 3-body parameter A*, one might 
worry that it may also depend on infinitely many other low-energy parameters associated 
with 4-body and higher n-hodj observables. If this is the case, the dilute BEC with large 
scattering length would have no universal properties at large na^. The most favorable 
possibility is that in the limit of large scattering length, the low-energy observables in the 

4- body and higher n-body sectors are all calculable to leading order in £/a in terms of a and 
A* only. If this is the case, then a dilute BEC with large na^ may have properties that are 
universal in the sense that they are determined by a and A^, only. 

This research was supported in part by DOE grant DE-FG02-91-ER4069 and by NSF 
grant PHY-9800964. 

Note added. — After this work was submitted for publication, we became aware of a 
related study by Bulgac [EO]. 
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